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Sediment transport is studied as a function of the grain to fluid density ratio using two phase numerical sim- 
ulations based on a discrete element method (DEM) for particles coupled to a continuum Reynolds averaged 
description of hydrodynamics. At a density ratio close to unity (typically under water), vertical velocities are 
so small that sediment transport occurs in a thin layer at the surface of the static bed, and is called bed load. 
Steady, or 'saturated' transport is reached when the fluid borne shear stress at the interface between the 
mobile grains and the static grains is reduced to its threshold value. The number of grains transported per 
unit surface is therefore limited by the flux of horizontal momentum towards the surface. However, the fluid 
velocity in the transport layer remains almost undisturbed so that the mean grain velocity scales with the 
shear velocity m*. At large density ratio (typically in air), the vertical velocities are large enough to make the 
transport layer wide and dilute. Sediment transport is then called saltation. In this case, particles are able 
to eject others when they collide with the granular bed, a process called splash. The number of grains trans- 
ported per unit surface is selected by the balance between erosion and deposition and saturation is reached 
when one grain is statistically replaced by exactly one grain after a collision, which has the consequence 
that the mean grain velocity remains independent of . The influence of the density ratio is systematically 
studied to reveal the transition between these two transport regimes. Based on the mechanisms identified in 
the steady case, we discuss the transient of saturation of sediment transport and in particular the saturation 
time and length. Finally, we investigate the exchange of particles between the mobile and static phases and 
we determine the exchange time of particles. 



I. INTRODUCTION 

After the_ pioneering works of RichardsorP, Rous^ 
and Vanoni'^', transport and dispersion of impurities sus- 
pended in turbulent flows, such as sand grains, dust, bub- 
bles or droplets, have received a renewed interest in the 
last decade, both from the fundamental point of view3fs! 
and for its applications to planetology-^, cloud physics® 
and geomorphology. In the later case, sediment may be 
entrained, transported and deposited by water flow or 
by wind. Then, gravity cannot be neglected as trans- 
port usually takes place in a turbulent boundary layer 
bounded by an erodible granular bed. Moreover, trans- 
ported particles are not passively advected by the flow: 
they induce a negative feedback, which eventually limits 
the erosion of the granular bed, leading to a steady state 
in which erosion and deposition balance each other. 

In such a homogeneous and steady situation, the fluid 
flow can be characterized by an unique quantity: the 
shear velocity m*, defined by the wall shear stress. The 
flux of sediments transported by the flow, called the sat- 
urated flux and noted Qsat, is an increasing function of 
Us, whatever the nature of the fluid. For aeolian trans- 
port, there has been great effort to obtain the rela- 
tion (?sat(M*) experiment all j^sHUl using both wind tun- 
nels and atmospheric flows in the field, numericalljl^SHSil 
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and theoreticalljff^'^lHlIl^ Similarly, several expressions 
for subaqueous bed-load have been proposed,^^ Most 
models are based on the same dynamical mechanisms and 
differ only by the approximation used to compute the 
particle trajectories^^ For this reason, experiments 
have been performed to determine the saltating motion 
of individual particles under wa.te^^^^ and in aiJ^^HUl 

Despite this wide literature, some fundamental as- 
pects of sediment transport are still only partly under- 
stood. For instance, the dynamical mechanisms limit- 
ing sediment transport, in particular the role of the bed 
disordeil^ and turbulent fluct uat ion^^SMU^ remain mat- 
ter of discussion. Also, derivations of transport laws 
have a strong empirical or semi-empirical basis, thus lack- 
ing more physics-related inputs. Here we investigate the 
properties of sediment transport using a novel numerical 
description of particle-laden flows, using two-phase nu- 
merical simulations based on a discrete element method 
(DEM) for particles coupled to a continuum Reynolds 
averaged description of hydrodynamics. In particular, 
we examine the transition from bed-load to saltation by 
studying the influence of the grain to fluid density ratio. 
A similar approach has recently been used to study the 
onset of aeolian saltatiorf22l. 

In section II, we introduce the equations of motion for 
the grains as well as the equations of hydrodynamics, em- 
phasizing the coupling between the two. Then, in section 
III, we detail the characteristics of saturated transport 
in the two limiting cases: bed load (water) and salta- 
tion (air) . In section IV, we propose an interpretation of 
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the simulations based on simple transport models. We 
then use these transport descriptions to derive and dis- 
cuss out-of-equilibrium transport and in particular the 
saturation length and time (section V). We contrast this 
time with the 'exchange time', which characterizes the 
diffusion of particles through the static/mobile interface. 
Finally, conclusions are outlined in the last section. 



II. TRANSPORT MODEL 
A. Key ideas 

We wish to model the transport of non-cohesive grains 
by a flow, under gravity. Although a continuum 'two 
phase' (grains and fluid) modeling^ is very appealing, it 
is problematic by several aspects, (i) It postulates that 
particles constitute an Eulerian phase, which means that 
the particles crossing an arbitrary control volume have 
almost the same velocity. In a homogeneous steady flow, 
an Eulerian approach immediately predicts that particles 
are transported along the direction parallel to the bed 
and to the flow - vertical velocities are ignored. However, 
at least for saltation, they are essential, (ii) Such a con- 
tinuum approach ignores the discrete and disordered na- 
ture of the granular phase. However, these properties are 
essential close to the transport threshold, below which no 
grain can be entrained. Such models incorrectly predict 
the threshold shear velocity and in particular its strong 
decrease with the grain Reynolds number. To avoid these 

issues, we use here a discrete element method for the 
particlePHsl. 

Resolving hydrodynamics around grains is technically 
feasible only if the size of the domain (the number of 
grains) and the time over which the simulation is run 
are very small. The idea introduced here is thus to use 
a continuum description of hydrodynamics, averaged at 
a scale larger than the grain size. This means that the 
feedback of the particles on the flow is treated in the 
mean field manner. 

This method allows one to perform very long numeri- 
cal simulations (typically 1000 ^/d/g), using a (quasi) 2D 
large spatial domain (typically 15000 spherical grains in 
a xyz box of respective dimensions 1000 dxldx 1000 d), 
while keeping the complexity of the granular phase. Pe- 
riodic boundary conditions are used in the x (fiow) di- 
rection. The domain height is large enough as to prevent 
any grain for reaching the top border. We will now de- 
tail the different ingredients of the model. To avoid the 
formation of ordered structures in the grain packing, we 
have used a slightly polydisperse sample (the maximum 
grain diameter is about 20% higher than the mean). For 
the sake of simplicity, we only give here the equations for 
the strictly monodisperse case (grains of diameter d). 



B. Forces on particles 

1. Equations of motion 

The grains have a spherical shape and are described by 
their position vector r, velocity u and angular velocity lj. 
A given grain labelled p inside a fluid obeys the equations 
of motion, 

q 

d^^rf (1) 
dt 2 ^ ^ ' 

where g is the gravity acceleration, / ~ md^/lO is the 
moment of inertia of a sphere, /P'* is the contact force 
with grain q, fp'''^ is the contact direction, and f^^[^ en- 
codes forces of hydrodynamical origin. 

2. Contact forces 

Following a standard approach for the modeling of con- 
tact forces in MD codes, see^"*^®^ and references therein, 
we consider the case where grains in contact are subject 
to (i) normal repulsion, (ii) tangential friction and (iii) 
energy dissipation. For simplicity, the normal repulsion 
is given by a spring-like elastic force - the results are 
independent of the spring stiffness, provided it is large 
enough. The tangential friction is modeled by a tangen- 
tial elastic force proportional to the relative tangential 
displacement between the grains. The moment of this 
force can induce particle rotation. Whenever the tangen- 
tial exceeds a given fraction of the normal force, defined 
by a microscopic friction coefficient, the contact 'slides' 
(Coulomb friction law). Finally, energy dissipation at the 
contact is ensured by adding a damping term to the force, 
proportional to the relative contact velocity. This term 
accounts for the restitution coefficient e, i.e. the ratio 
between grain velocities after and before a collision. 

3. Hydrodynamic forces 

For simplicity we assume that the net hydrodynamical 
force (/^ui(j) acting on a grain p due to the presence of 
the fluid is dominated by the drag and Archimedes forces, 
/lag and fl^.^^^, respectively: 

/fluid ^ /drag + /Arch (2) 

The lift force, lubrication forces and the corrections to 
the drag force (Basset, added-mass, Magnus, etc.) are 
neglected. Notice that an important consequence of ne- 
glecting lubrication forces is the constancy of the grain's 
restitution coefficient, which in general is function of the 
grain's Stokes number 



3 



a. Drag force: We hypothesize here that the drag 
force exerted by a homogeneous fluid on a moving grain 
only depends on the difi^erence between the grain velocity 
iPlXjz) and the fluid velocity u(z) at grain's height z. 
Introducing the particle Reynolds number based on 
this fluid-particle velocity difference i?„ = \u — iP\d/v, 
the drag force can be written under the form 



•J drag 



-pf(fCd{Ru)\u - vF\{u - vP 



(3) 



where Cd{Ru) is the drag coefflcient and p/ is the density 
of the fluid. We use the following convenient phenomeno- 
logical approximatiorPsl 



Cd{Ru) — 



(4) 



where ~ 0.5, is the drag coefflcient of the grain in 
the turbulent limit ^ oo), and R^ ~ 24 is the tran- 
sitional particle Reynolds number above which the drag 
coefficient becomes almost constant. 

b. Archimedes force: This force results from the 
stress which would have been exerted on the grain, if 
the grain had been a fluid. Thus, 
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(5) 



where ^(P is the e;ram volume and a 



„ . T^- is 

the undisturbed fluid stress tensor (written in terms of 
the pressure and the shear stress tensor tJj). In flrst 
approximation, the stress is evaluated at the center of 
the grain. 

In the inner region of the boundary layer, where 
we assume that most of the transport takes place, the 
Archimedes force can be well approximated by: 



J A 



Arch - 



(6) 



where Cx and Cz are the horizontal and vertical directions, 
respectively. 



C. Hydrodynamics 

In the presence of particles occupying a volume fraction 
(/), the hydrodynamics is described by the two-phase flow 
Reynolds averaged Navier-Stokes equations: 

Pf{l - 4>)l}tu, = -d,pf + pf{l - cj))g, + a,r^. - (7) 



where DtUi = dtUi+UjdjUi denote the fluid inertia, t/j is 
the total shear stress tensor resulting both from viscous 
diffusion of momentum (viscous stress) and transport of 
momentum by turbulent fluctuations (Reynolds stress). 
F is the body force exerted by the grains on the fluid. 
It reflects the velocity fluctuations induced by a moving 
grain. As we focus in this paper on steady homogeneous 
sediment transport, we hypothesize that the influence of 



a given grain remains localized in a thin horizontal region 
and that the typical horizontal distance over which the 
flow is disturbed is comparable to the distance between 
moving grains. F(z) can then be obtained by averaging 
the hydrodynamical force /g^jj acting on all the grains 
moving around altitude z, in a horizontal layer of area A 
and thickness dz: 



F{z) 




(8) 



Here, we take for A the total horizontal extent of the do- 
main (i.e. 1000c? X Id). Notice that the number of grains 
given by the condition p £ {z; z-l-dz} is also a differential 
quantity. The symbols (.) denote ensemble averaging. In 
order to gain statistics, we make use of the steady charac- 
ter of the studied situation, and also use time averaging. 
For simplicity, wc note r-' = t^^ the fluid shear stress, 
and u = Ux for the fluid horizontal velocity. After sub- 
stituting the drag and Archimedes forces explicitly, F[z) 
can be rewritten as, 



F{z)^4>iF^, 



0(a 



- dzp- 



(9) 



where the grain's volume fraction (/) is defined as 

^3 (10) 



^ ' A<\z ^ 6 

p€{2;z+dz} 



and 



^drag)(^) 



E /lag)/ E 



I p€{2:z+dz} 



pG{2;z+dz} 



(11) 

is the average drag force acting on grains at height z per 
unit grain's volume. 

In the inner region of the turbulent boundary layer, 
both the fluid inertia and the horizontal stress gradi- 
ents can be neglected, and the vertical component of the 
Reynolds equation becomes 

O--az/-p/(l-0).g-i^z (12a) 

or dzP^ = -pfQ - YZ^^ (^drag.z) (12b) 

Under the assumption of steady and homogeneous sedi- 
ment transport, the contribution of grains' vertical drag 
to the momentum balance is negligible and the vertical 
balance reduces to the hydrostatic pressure: 

9z/ = -Psg (13) 
and thus the Archimedes force ([6]) simplifies to 

/Arch = {dzT^Cx - Pfg) (14) 

which reduces to the buoyancy force for a static fluid 
(9zr^ = 0). 
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Furthermore, after neglecting inertia and horizon- 
tal stress gradients, the horizontal component of the 
Reynolds equation becomes 



(15a) 



which, with the use of the x-component of Eq. [9j can be 
rewritten as 



1 



drag, X / 



(16) 



The prefactor (1 — (j>)~^ accounts for the increasing role 
of the Archimedes force at the bed, where the volume 
fraction reaches its maximum value 4>b- Notice however 
that this increase is balanced by the decrease of the aver- 
age drag force as the fluid velocity approach zero inside 
the bed. 

For simplicity, in what follows we retain the grain's 
horizontal feedback term F^, including Archimedes' 



force, and integrate (15a I as 



{z) — Pfui — T^{z) 



(17) 



where we have introduced the shear velocity w*, defined 
by the undisturbed (grain free) wall shear stress, and the 
grain borne shear stress defined by 



F,iz')dz' 

A, 

(^drag, x) dz' 



1 



(18a) 
(18b) 



where we have replaced dzT^ by (16). In the integration 



bounds, CXI means that the integral includes all moving 
grains. Here we are assuming that grain's motion takes 
place well inside the boundary layer, such that at the 
upper limit of the transport layer, where F^ = 0, the 
horizontal shear stress is constant and equals Pfu^. 

In order to relate the fluid borne shear stress to the av- 
erage fluid velocity field, we adopt a Prandtl-like turbu- 
lent closur^^. Introducing the turbulent mixing length 
£, we write 



rf - 



(19) 



v is the viscosity (a constant independent of the vol- 
ume fraction). As for the mixing length £, we know it 
should vanish below some critical Reynolds number 
and should be equal to the distance to the surface z, far 
above the transport layer. A common phenomenologi- 
cal approach is to express the turbulent mixing length 
as a function of the Reynolds number and z (see below 
Eq. 21 ). However, this involves the definition of an inter- 
phase between the static and mobile zones, below which 
£ must vanish. To avoid the need of such a somewhat ar- 
bitrary definition, we propose instead a differential equa- 
tion 



dj = K 



1 — exp — 




(20) 



where n ~ 0.4 is von Karman's constant and the dimen- 
sionless parameter i?c = 7 is determined from an indi- 
rect comparison to measurements (see below). The ratio 
u£/v is the local Reynolds number based on the mixing 
length. Note that a function other than the exponential 
can in principle be used, provided it has the same behav- 
ior in and — oo, although the present choice provides a 
quantitative agreement with standard data. As discussed 
below, this formulation allows us to define £ both inside 
and above the static granular bed. Note that, although 
it appears here as a side technical aspect of the present 
study on sediment transport, Eq.|20]is in fact an original 
piece of work, interesting in itself, and which can be used 
in other contexts. 



D. Solutions for special cases and validation 

Before presenting the steady state properties of the 
transport model, let us first provide some technical de- 
tails of the integration of the model and compare some of 
the results to what is known for two simple cases: the law 
of the wall for smooth and rough surfaces. We will then 
briefly discuss the hydrodynamics for stationary trans- 
port. 



1. Smooth surface 



In the absence of grains, — Pful and the flow ve- 



locity u{z) is given by (19), with the boundary condi- 
tion m(0) ~ 0, and some mixing length phenomenological 
equation. A popular choice for the flow over smo oth su r- 
faces is arguably the one suggested by van DriestPl^ 



1 — exp 



1 



i?vD V 



(21) 



which reproduces well classical experimental results, in 
particular the transition from the viscous sub-layer where 
the velocity profile is linear u(z)/u* = zu^jv^ to the tur- 
bulent region, where the velocity profile is logarithmic 
u{z)/u^, — K~^logz/zo. In this smooth case, the hydro- 
dynamic roughness length zq scales on the viscous length: 

— 0.1j^/u*. Quantitative agreement with this scaling 
law is obtained with _RvD — 26. 

In contrast to van Driest's mixing length, which can 
be directly computed for any distance z to the surface, 
the use of the differential equation (20) necessitates to 
start from a boundary condition at z = 0. Expanding 
dz£ around u£ = Q and substituting the viscous velocity 
profile u = zul/v, we can get the asymptotic behavior of 
£{z) close to the surface and then use it as the boundary 
condition for z — 0. The parameter i?c = 7 is then fitted 
to achieve the best comparison of the velocity profile ob- 
tained by integration of Eq. [20| to the one obtained using 
the van Driest's expression (21 ). Figure [l] shows that this 
comparison is quantitative. 
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FIG. 1. Comparison of the normalized vertical velocity profile 
using the differential mixing length equation proposed here 
( 20 1 with Rc — 7 (long-dashed line) and the standard van 
Driest equation \2\\ with _RvD = 26 (solid line). The two 
asymptotic behaviors are also shown: the linear velocity pro- 
file at the viscous sub-layer (dashed line) and the logarithmic 
one far from the surface (dotted line). 



2. Static bed 

The real advantage of Eq.[20]is evident in the presence 
of a granular bed. In this case, the interface between the 
bed and the fluid is not well defined, making the use of 
Eq. [21] difficult. However, as shown below, Eq. [20] can 
be integrated from an arbitrary position inside the bed, 
where the asymptotic solution [z — > —00) holds, so that 
the result of integration is independent of that position. 

a. Asymptotic velocity profile: Let us consider a 
static {uP = 0) and homogeneous bed with an average 
volume fraction <j)\,. As the flow velocity inside the bed 
is very small, the average drag force per grain's volume 
can be approximated by the viscous drag (Fdrag, a;) = 
ISpfVu/cf. Besides, the mixing length is also very small 
and the fluid borne shear stress is dom inated by the vis- 
cous diffusion = pjudzU. From (16 1 it follows: 



a u 



(t^b 18 

1 - 06 d2 



u{z) 



(22) 



The solution of this equation is an exponential decay of 
the velocity inside the bed with a characteristic 'penetra- 
tion' length 



Afc = d^{l - (/.6)/(1806) 



(23) 



In our case, we use spherical bed with </ib ~ 0.6 and thus 
A ~ 0.2d. Integration of the velocity can be started at 



FIG. 2. Dimensionless roughness length for increasing ratios 
of grain diameter to viscous length du,/v (symbols). The 
prediction for smooth surfaces zo ~ Q.\v/u, (full line) is also 
shown for comparison. 



any arbitrary position well inside the bed (typically 4- 
5 grains depth is enough) . The starting value is found 
iteratively to ensure that the first integration steps follow 
the asymptotic behavior 



u{z) ^ Us exp ((z - Zs)/Xb)- 



(24) 



b. Asymptotic mixing length profile: As both the ve- 
locity and the mixing length are vanishingly small inside 
the bed, ( 20 ) expands as 



(25) 



which, after substituting the exponential velocity profile, 
integrates as 



'-u{z) 



(26) 



The proportionality between i and u is then used to set 
the initial value for the mixing length well inside the bed. 

c. Solution for a static bed: From velocity profiles, 
we see two characteristic length scales: the penetration 
length A ~ 0.2d determined by the bed volume fraction 
and the viscous length v/u^ in the viscous sub-layer. For 
the numerical integration of the velocity we use a fraction 
of the smallest scale of the two as the mesh size dz. The 
mesh is linear inside and just above the bed. Few grains 
above it, we proceed integration by exponential vertical 
steps, to adjust better to the logarithmic velocity region. 
After integration, the flow velocity can be computed at 
any position by either linear or logarithmically interpo- 
lation. 
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Figure [2] shows the hydrodynamical roughness length, 
obtained from the logarithmic velocity profile far from 
the bed surface, for different ratios of the grain diame- 
ter d to the viscous length v/u^. There are two distinct 
regimes: (i) when the viscous length is on the order or 
higher than the grain diameter, the granular surface be- 
haves as hydrodynamically smooth, with the roughness 
length proportional to the viscous length zq oc u/u^,; (ii) 
when the grain diameter is much larger than the vis- 
cous length, there is turbulent dissipation at the surface 
and the bed becomes hydrodynamically rough with the 
roughness length seeming to converge to the value around 
d/5. Note that this value of the roughness is larger than 
what is classically reported from experiments {zq ~ d/30 
iiPl zo ^ d/24: iiPl zo ^ d/10 in'^ we don't expect 
quantitative agreement working with a two-dimensional 
model, and model parameters could be more refined for 
this purpose. 



3. Mobile bed 
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When grains of the bed move, their motion is confined 
to the surface and the fluid velocity profile is still expo- 
nential inside the bed. Therefore, we can use the same 
procedure to integrate the equations as for a static bed 
above. As an example, Fig. [3] shows the hydrodynamic 
profiles for the steady state of a transport simulation with 
density ratio Pp/ Pf = 2. In what follows, as a convention 
for all vertical profiles, the reference height z = is by 
definition the height at which the volume fraction equals 
half the bed one, 0(0) = 0h/2. Notice this definition is 
not related to grain motion, and finite grains velocity are 
effectively observed for negative values of z. 

Once the grains are able to move, a transport layer 
quickly develops, accompanied by a strong negative feed- 
back on the flow, and the subsequent reduction of the 
fluid borne shear stress towards zero inside the bed 
(Fig. [3]). As expected, close to the bed (typically z < d) 
the viscous shear r/ = pfvdzU is dominant as the mixing 
length is very small (about O.Old). Notice that all profiles 
are smooth and continuous across the bed surface: (i) in 
the static bed {(j)/(pb — 1) there is no turbulence and the 
balance of viscous drag and viscous momentum dissipa- 
tion leads to the exponential decay of the fiow velocity; 
(ii) at bed surface {4'/4'b ^ 1) grains within the transport 
layer are accelerated until the momentum dissipated due 
to collisions balance the momentum extracted to the fiow; 
and (iii) at the undisturbed boundary layer {(j)/(j)b — >■ 0) 
the fiuid shear stress is constant and the flow momen- 
tum is dissipated through viscosity and mostly turbulent 
mixing (Fig. [s]). 

E. Dimensionless numbers for transport 

We now discuss how to make the model equations di- 
mensionless. Gravity gives the relevant scale for forces. 



FIG. 3. Vertical profiles of the rescaled volume fraction 
(l>/(l>b, flow velocity = u/u*, mixing length — £/d, 
fluid borne shear stress — /pfU^, viscous shear stress 
r/^ = vdzu/u^ and turbulent shear stress r/^ = {Idzuf' /u1 
(by definition t^'^ = r/^ -I- t/'^). The reference height 2 = 
is set at the altitude such that — 4'b/2. 



TABLE I. Units used in the model, expressed in terms of the 
grain density (pp), fluid density (p/), gravity [g) and mean 
grain diameter (d) 



General 


length I 


d 


acceleration 


9 


time t 




velocity v 


Vgd 


Particles 


angular velocity uj 




mass m 


IPvd' 


moment of inertia / 


md? 


force / 


mg 


contact stiffness k 


mg/d 


damping constant 7 


ra\J gjd 


Fluid 


shear stress r 


{Pv - Pf)9d 



More precisely, it only appears in the grain equation 
of motion under the form of a buoyancy-free gravity 

^1 — ^ j g. The choice of the typical length scale is less 

obvious. On the one hand, the contact forces and the 
trapping of particles at the surface of the bed do not 
depend on the fluid properties: the grain diameter d is 
thus the relevant length scale for the static grains. On the 
other hand, one can build a drag length from hydrody- 



namics, which is the length needed to accelerate a grain 
to the fluid velocity. This inertial length is proportional 
to — c? and is the relevant length scale for the mobile 

grains. This means that the density ratio Pp/pf cannot 
be eliminated and is a true dimensionless parameter of 
the problem. We shall see below that this density ratio is 
the parameter controlling the transition from bed load to 
saltation. We have chosen d as a reference length scale, 
and Table |T] summarizes all the parameters used in our 
code. The second control parameter is the shear velocity 
M* imposed far from the bed, or equivalently the shear 
stress Pfyj. Its dimensionless counterpart is the Shields 
numbeiES, defined by 



e = 



Pfi 



{Pp - Pf )9d 



(27) 



which encodes the strength of the flow. Making 
the viscosity non-dimensional, we obtain a grain-based 
Reynolds number 



R. = 




-l]gd 



(28) 



Physically, it determines the hydrodynamic regime at the 
scale of the grain. The figures presented in this paper 
are obtained using the same particle Reynolds number 
Re = 10. This value is sufficiently large to ensure that 
the grain diameter is about one order of magnitude larger 
than the viscous length. As a consequence, the flow is 
fully turbulent at about one or two grain layers above 
the static grains. 

Dynamics at the scale of the contact between grains is 
controlled by other numbers: the (here constant) resti- 
tution coefficient e = 0.9, the frictio n coefficient p = 0.5 
and the contact duration — tt/ y/2k/m — (7/m)^ (in 
terms of the contact stiffness k ~ 5000 mg/d and the 
damping constant 7 = ^,.1 m^fgjd, se^^. We have 
checked that the values given to these parameters do not 
qualitatively change the results. 



III. SATURATED TRANSPORT 
A. Qualitative results 

Transport equations are integrated until a statistically 
steady homogeneous state is reached. Since we are pri- 
marily interested in the transition from bed load to salta- 
tion, we have varied the Shields number within the range 
= 0.003-0.5 (a range which contains the threshold 
0d, see below) and the density ratio within the range 
p^lpj = 2-2000. 

Once transport has reached its saturated state, the 
general picture is as follows: at small density ratios 
Pp/pf — 2, which is the typical value under water, the 
transport is confined at the surface, within a few grain di- 
ameters. The dense and thin transport layer is character- 
istic of the bed load regime. On the contrary, at a large 
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FIG. 4. Transport profiles: volume flux density 19(2;) (left) 
and volume fraction (l>{z) (right) for water (solid lines) and 
air (dashed lines). 



density ratio Pp/pf — 2000, which is typical of aeolian 
situation, the transport layer becomes wide and dilute, 
extending over several tens of grain diameters (Fig. |4|. 
This is typical of the saltation regime. Within the very 
same numerical model, we are thus able to reproduce the 
basic characteristics of transport at both limits. 



B. Saturated flux 

Steady and homogeneous sediment transport is basi- 
cally quantified by the volumetric saturated fiux (/gat, i-e. 
the volume of the particles (at the bed density) crossing 
a vertical surface of unit transverse size per unit time. It 
has the dimension of a squared length per unit time. In 
the simulations, we compute it as 



9sat 
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p 



(29) 



A key issue is the dependence of (?sat on the shear velocity 
or, equivalently, on the Shields number O. In order to 
highlight this dependence, figure [5] shows the saturated 
fiux rescaled by 9 in both cases (wate r and air) . In ag ree- 
ment with experimental observation ji^ l 32H36 | 39 | 76| t79l .^g 

find that ggat scales asymptotically as O (or u^) for 
saltation, while gsat scales as O^/^ (or uj) underwater 
(Fig. [5]). Most models of aeolian transport miss the in- 
fiuence of the negative feedback of transport on the fiow. 
Therefore, they do not give the correct scaling, predicting 
Qsat oc . We demonstrate below, in the same numerical 
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FIG. 5. Rescaled saturated flux Qsat/Q^/ {pp/ pf — i)gd'^ ver- 
sus the rescaled shear velocity y/Q for water (a) and air (b). 
For air the saturated flux scales asymptotically as while 
for water it follows B'^''^. Full lines are the predictions of the 
simplified models for bed load (Eq. |40[ ) and saltation (Eq. |46[ ), 
given in the text. 



model, a fundamental difference between the two trans- 
port regimes, which correspond to different underlying 
dynamical mechanisms. 

Figure [5] reveals the existence of a threshold shear ve- 
locity below which the flux vanishes. More precisely, we 
define the dynamical threshold Shield number from 
the extrapolation of the saturated flux curve to 0, which 
gives in our case 8^ ~ 0.12 for water {pp/pf = 2) and 
Qd — 0.004 for air {pp/pf = 2000), respectively. These 
values are consistent with experimental ones within a 
factor of 2. Once again, a refined tuning of these val- 
ues could be achieved by adjusting the model parameter 
(e.g. Rc) and performing 3D simulations. Figure [6] shows 
the dependence of the threshold with the density ra- 
tio pp/pf. It is usually assumed that the Shields num- 
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FIG. 6. Dynamical threshold Shield number Qd as a function 
of the density ratio pp/pf- 



ber compares directly the horizontal force exerted on a 
surface grain to its weight, in which case the threshold 
Shields number could be interpreted as an effective fric- 
tion coefficient, within a numerical factor. If this was 
true, Qd would be a constant, independent of Pp/pf- 
However, one observes that Qd decreases rapidly with 
the density ratio. 



C. Transport layer 

Figure [7] presents the vertical profiles of the flux den- 
sity, i.e. the flux per unit height (such that qs^t = 
J ■d{z)dz) for different shear velocities. It shows that bed 
load and saltation mainly differ by the vertical charac- 
teristics of the transport layer. At small density ratios 
the motion of grains is confined within a thin layer of few 
grain diameters (Fig. [7^) . Most of the bed load occurs at 
about one grain diameter above the static bed and the 
flux density profile decays symmetrically on both sides 
of this maximum. By contrast, for large density ratios, 
grains experience much higher trajectories and the trans- 
port layer is much wider. Figure [7)d shows that the flux 
density still presents a maximum close to the static bed 
but decreases exponentially with height. 

These qualitative observations can be formalized by 
deflning a characteristic transport layer thickness A from 
the flux density proflle i}{z) as: 



7„°°(z-z)2^(z)dz^ 



1/2 
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(30) 



where , 



z'd{z)dz gives the altitude of the trans- 
port layer centre. If the flux profile decreases exponen- 
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FIG. 7. Vertical profiles of the sediment flux density ^{z) for 
different values of the shear velocity ratio ^/O/S^, in water 
(a) and air (b). Insets: characteristic transport layer thickness 
A as function of the shear velocity. 



tially, A is the characteristic distance over which this de- 
crease takes place. The variations of A with the shear 
velocity are presented in the insets of figure [Tj For un- 
derwater bed load, the size of the transport layer is about 
one grain diameter, gently increasing with the shear ve- 
locity from A ~ d/2 to A ~ d. For aeolian saltation the 
transport layer is indeed wider, with a characteristic size 
A ~ 50d roughly independent of the shear velocity. 

Figure |8] shows the dependence of the transport layer 
thickness A with the density ratio. At large density ratios, 
A is observed to scale with ^d, i.e. proportional to the 
drag length, which emerges when the motion of the grains 
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FIG. 8. Characteristic transport layer thickness A (•) as func- 
tion of the density ratio for B = 20d- At small density ratios 
it is limited by the grain size (dashed line), while for large 
ones it scales as pp/pfd (solid line). The averaged vertical 
energy per grain {{w^)^)/g (o) is also shown to illustrate the 
dynamical origin of A (see text). 



is dominated by the balance between inertia and hydro- 
dynamical drag. This length is thus expected to control 
the characteristic hop height and hop length, which natu- 
rally leads to wider transport layers for lighter fluids. In 
sub-aqueous conditions (density ratios of order unity), 
the transport layer thickness is limited by the grain size 
A ~ d, which is the characteristic length sc ale for contact 
forces and geometrical trapping of particle^snilll_ 



The hop height can be estimated from the particle ver- 
tical velocity w'^ using the ballistic approximation, ne- 
glecting the vertical component of the drag force. Under 
this hypothesis, one expects the hop height to increase 
like {w^Y Ig. Figure [s] shows the dependence of the av- 
erage squared vertical velocity {{w^Y) on the density ra- 
tio. One observes that the transport layer thickness A 
is effectively determined by the hop length {{w'^)'^) / g for 
Ppl Pf ^ 10- Below this cross-over value, the transport 
layer thickness is given by the grain diameter d, as tra- 
jectories are almost horizontal. The transition from bed 
load to saltation therefore takes place when the vertical 
velocities of the particles are sufficiently large for these 
particles to escape the traps formed by the grains on the 
static bed. Formally, the criterion of this transition can 
then be written as ((w^)^) ~ gd. 
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D. Grain feedback on the flow and the conditions for 
saturation 

Another difference between bed load and saltation is 
how the grain's feedback on the flow is distributed within 
the steady state transport layer. Figure [9] presents the 
vertical profiles of the fluid shear stress, rescaled by the 
dynamical threshold = Qd{Pp/pf — ^)gd (as defined 
by the saturated flux), for different shear velocities. For 
bed load (Fig. [9^), the different profiles of the fluid shear 
stress seems to converge to the threshold value very close 
to the surface (z = 0). As the dynamical threshold repre- 
sents the limit shear stress for the fluid to sustain trans- 
port, the grain motion is directly driven by the fluid and 
thus controlled by the excess shear stress Pful — r^. In 
this transport layer, the fluid momentum decays over few 
grain sizes, in agreement with the vertical extension of 
the transport layer. In contrast, the fluid shear stress is 
below the threshold in the bed {z < 0) but some (weak) 
transport still occurs there, which is sustained not by 
the fluid itself but by the momentum transferred to the 
surface by grain collisions. 

This general picture is still valid for saltation (Fig.|9]3), 
however now the dynamical threshold is reached much 
farther from the surface (at z ~ 10 c?) which implies that 
the kinetic energy of impacting grains is large enough as 
to sustain the transport below this height. Above it, the 
transport is driven by the fluid and most of its momen- 
tum is dissipated in a much larger layer (comprising tens 
of grain diameters) again in agreement with the size of 
the saltation layer. Notice that although this surface sub- 
layer below 10 d contains most of the grains (see Fig. |4]) 
it still represents a small fraction of the overall transport 
layer, whose characteristic thickness is A ~ 50 d. 

An important consequence of this distinction in the 
vertical structure of the grain's feedback is that although 
for bed load transport is equilibrated when the fluid shear 
stress reaches its dynamical threshold below the trans- 
port layer, this condition is not enough for saltation to 
equilibrate. For saltation there is a sub-layer where trans- 
port is not directly driven by the fluid and thus its equi- 
libration is not dictated by the threshold. There, the 
properties of grain's collisions become relevant and the 
equilibrium is described by the conservation of the num- 
ber of saltating grains i.e. when the number of grains 
entering the flow exactly balance those grains trapped 
by the bed. 




FIG. 9. Vertical profiles of the fluid borne shear stress t^(z) 
rescaled by the dynamic threshold Td — Qd{pp/pf — ^)gd for 
different values of the shear velocity ratio ^/O/Od, in water 
(a) and air (b). 



IV. INTERPRETATION 

A. A simple transport model for bedload 

We propose here a simple model of bed load inspired 
from Bagnold's original idea^^. We hypothesize that 
moving grains are confined in a thin layer with thickness 
on the order of d. As the average particle vertical ve- 
locity is very small, grain hop heights are typically much 



smaller than d (Fig. [S]), which means that the vertical 
motion of the grains can effectively be neglected. The 
saturated fiux can then be decomposed as the product of 
the number n of transported grains per unit area by the 
mean grain horizontal velocity vP: 
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(31) 
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In the numerical simulations, we compute n and as 



2 

1 



(32) 
(33) 



Saturation is reached when the fluid shear stress equals 
the transport threshold at the surface of the static bed 
i.e. when = pju\ ~ r^, with = Qa{PplPf - ^)gd 
(Fig. [9^). As consequence, the number of transported 
particles per unit area is solely determined by the excess 
shear stress: 



p "p 



Notice that these definitions are consistent with the def- 



fd 



lPd<P 



inition of g,at If all grains were moving at the same YinaWy, the saturated flux reads: 



velocity, then n and u'^ would indeed be respectively the 
density of moving grains and their velocity. 

We can then write the grain born shear stress as pro- 
portional to the moving grain density n and to the drag 
force acting on a grain moving at the average velocity vP 
due to a flow at the velocity u: 

rP^nfd with fd = lcTpf{u-uPfd\ (34) 

o 

For the sake of the argument, we neglect Archimedes con- 
tribution as well as the dependence of the drag coefficient 
on the particle Reynolds number. A key assumption is 
that grains are in a steady motion, which means that 
the drag force fd balances a resistive force due granular 
friction, collisions with the bed, etc. These different dissi- 
pative mechanisms can be modeled as an overall effective 
friction force characterized by a friction coefficient /i^: 



fd = ^PdiPp - Pf )9d^ 



(35) 



We can furthermore express the fluid velocity Ud at the 
transport threshold by assuming that the hydrodynamic 
drag exerted on a static grain (u^ = 0) has to overcome 
a static friction, characterized by a coefficient fig'. 
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(40) 



Inserting the expression (36) oi Ud, one gets the scaling 
law for the flux at large O: 
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(41) 



B. A simple transport model for saltation 

We now proceed in a similar manner for the aeolian 
saltation regime, following ideas initially proposed by 
Owen (1964) and Ungar & Haff (1987). In this regime, 
the motion of the grains is not conflned to a thin layer 
at the surface of the bed. We consider an average grain 
trajectory, in which the particle takes off from the bed 
with the horizontal velocity u^, and comes back to it with 
a velocity u'^, after a hop of length a. Some momentum 
is extracted from the wind flow by the grains to perform 
their jumps, so that the particle shear stress writes 



Ud 
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(36) 



Combining the above equations it follows that the veloc- 
ity difference between the grain and the flow is constant: 



(37) 



We now assume that the transported grains do not dis- 
turb the flow. Then, the flow velocity around grains 
u must be proportional to the shear velocity, so that 
u/ud — yjQ/Qd (see Fig. 10 1). One therefore deduces: 




(38) 



This predicts that the grain velocity does not vanish 
at the threshold, if friction is lowered during motion 
(/id < ^s)- The velocity at threshold iid(l — ^ fid/Ps) 
can be interpreted as the velocity needed by a grain to 
be extracted from the bed and entrained by the flow. 
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(42) 



Now we use again the decomposition of the saturated 
flux as the product of the grain density n and the grain 
velocity (Eq. 31). Saturated transport corresponds to 
the balance = P/""* ~ Td, so that n still has the same 
form as in the bed-load case: 



(Pp - Pf)9d 

fd 



(0-0d): 



(43) 



but with a different effective drag force fd, not related to 
friction anymore but to grain velocities. As the grain hop 
length can be related to the grain velocity as a oc u^w^ /g 
(ballistic approximation), we can effectively write 



fd OC -d' 
6 



Pp9 



(44) 



Now, for saltation, steady transport also implies that 
the number of grains expelled from the bed into the flow 
exactly balances those trapped by the bed, i.e. a replace- 
ment capacity equal to one. Due to the grain feedback on 
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the flow, in contrast with bed load, grains in the trans- 
port layer feel a flow independent of the wind strength 
(see Fig. ^ and Fig. [To|3). Thus, new moving grains 
come only from high energy bed collisions. Since the 
number of ejected grains is a function of the impact en- 
ergy (or equivalently, of the impact velocity), the mean 
grain velocity must be constant, independent of the 
shear velocity, scaling with u^: 



CX Ud- 



(45) 



From this argument, it follows that all particle surface 
velocities (uj^, Up lu^) also scale with Ud, so that fd is 
also a constant. Finally, the scaling law followed by the 
saturated flux becomes, 



9sat oc (1 - pf/pp)udd (6 - Qd) ■ 



(46) 



C. Comparison with simulations 



The above simplifying models suggest a few simple 
tests to investigate the dynamical mechanisms in the 
DEM simulation, (i) Is saturation of transport due (or 
not) to the negative feedback of moving grains on the 
fluid? (ii) Do we recover the linear relation between the 
grain density n and the excess Shield number — Qd, 
whatever the transport regime? (iii) Does the mean grain 
velocity vP depend (or not) on the shear velocity? 



1. Grain feedback on the flow 

The information of the feedback of moving grains on 
the fluid flow is formally encoded in the flow roughness 
length zo- However, it can flrst be qualitatively under- 
stood from the shape of the fluid shear stress and flow 
velocity proflles inside the transport layer (Fig. |9] and 
Fig. 10 respectively). For bed load, as shown in the in- 
set of Fig. [TO^ , the flow velocity at the point of maximum 
transport (at about Id), increases with the shear veloc- 
ity. This indicates that the flow is barely disturbed in 
the transport layer. As shown in Fig. |9^, most of the 
feedback is actually conflned very close to the static bed 
and thus, most moving grains feel the undisturbed flow 
that increases with u*. In contrast, for aeolian transport 
(see Fig. [9]d and Fig. \lO\3) the flow velocity is strongly 
affected by the motion of grains as it becomes almost 
independent of the wind in the region z < lOd. 

The data of the hydrodynamical roughness length zq 
show a similar picture (Fig. 11). In the saltation regime 



the roughness length increases with the shear velocity as 
a result of grain feedback, which can be modeled from 
the existence of a focal point where u = Ufatz — Hf 
independently of (Fig. \lO\p) , and above which the flow 
velocity recovers its log proflle u ~ u^,/k \n{z / zq^^^^^^^. 
This gives 
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FIG. 10. Flow velocity vertical profiles at different shear ve- 
locity ratios ^/OjOd for water (a) and air (b). Insets: velocity 
at the point of maximum transport {z ~ Id), as a function 
of the rescaled shear velocity. The dashed line in the upper 
inset corresponds to the fit u oc u,. 



This expression reproduces well the increase of Zq for 
stronger winds, when the density ratio pp/pf is large 



enough (Fig. 111. Typically below Pp/pf — 10, Eq. 



47 



does not reproduce the data anymore. This is consistent 
with the absence of a focal point in the bed load regime 
(Fig.~ 



lengt 



10 1). Also, in the small Pp/pf limit, the roughness 



1 remains very small (substantially smaller than d). 



2. Number of transported grains and average grain 
velocity 



Zq ~ Hf exp {—KUf/u^,). 



From expressions ( 32 ) and ( 33 1 , we can compute the 



(47) number of transported grains per unit area and the mean 
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FIG. 11. Rescaled hydrodynamical roughness length as func- 
tion of the shear velocity for different density ratios. Solid 
lines are the predictions based on the focal point assumption 
(Eq.|47|. 



grain horizontal velocity as a function of the shear veloc- 
ity of the flow. Figure [T2] shows a linear relation between 
n and Q — Qd for both bed load and saltation. This is 
consistent with the predictions of the above models. In- 
terestingly, the friction coefficient /i^, defined from the 
proportionality factor (see Eq. 39), has the same value 
~ 1 in both cases. This suggests that dissipation due 
to collisions of the moving grains with the bed plays the 
same role in both transport regimes. 

The dependence of the mean grain velocity is also 
fully consistent with the picture emerging from the simple 
models. As shows in Fig. [TS} increases linearly with 
y/Q/Qd for bed load (Eq. 38 ) while it remains roughly 
constant for aeolian saltation (Eq. 45 1. Interestingly, the 
different curves shown in Fig. [13] cross at ~ Qd- In 
other words, the grain velocity at the transport thresh- 
old scale on \/gd{pp/ pf — 1), with a prefactor slightly 
smaller than unity, whatever the transport regime. This 
common behavior between bed-load and saltation results 
from the fact that the negative feedback of transport on 
the flow disappears at the threshold, as n vanishes. 

Fitting the grain density and the mean grain velocity 
to the simple model of bed load, one can extract the ef- 
fective friction coefficients pd and Ps- The static friction 
coefficient ps turns out to be ~ 6 times larger than the 
dynamical friction coefficient fid- This means that the 
motion is lubricated by the fluid once the grains are en- 
trained. Therefore, the grain velocity at the threshold 
remains finite (but n vanishes). 



FIG. 12. Linear relation between the rescaled number of 
transported grains per unit area and the rescaled excess of 
shear stress for different density ratios. 
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FIG. 13. Rescaled mean grain velocity as function of the 
rescaled the shear velocity for different density ratios. Full 
lines show the analytical prediction given in the text for the 
two limiting cases: water (Eq. 38 1 and air (Eq. 451. 



V. SATURATION TRANSIENT 

Beyond the properties of steady and homogeneous 
transport, we address in this section the time and length 
scales involved in the relaxation of the sediment flux to- 
ward its saturated value, which are relevant in the con- 
text of pattern formation^* " We further emphasize the 
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difference between the saturation time and the exchange 
time. 



A. Saturation length and time 

Whatever the transport regime, the saturation tran- 
sient is controlled by two mechanisms. On the one hand, 
bed erosion or deposition must take place to adapt the 
number of transported grains to the flow velocity. On 
the other hand, grains must be accelerated by the flow 
to their asymptotic velocity. 

We have addressed the case of saltation in a s eries of 
articles, starting from a controversy between u^SHSHH 
and resolving i1p2l. In summary, the horizontal accelera- 
tion of a grain entrained by the wind is governed by the 
equation of motion: 



duP 



4 ppd 
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(48) 



Contrarily to bed load, in the saltation regime, dissipa- 
tion only takes place during collisions and not through 
a permanent friction on the static bed. The only length 
scale in this equation is the so-called drag length y^d. As 
a consequence, the relaxation of the particle velocity to 
the fluid velocity occurs over a length which varies as 



Pp 
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independently of the wind speed, with a proportionality 
factor that depends on the restitution coefficient Ex- 
cept in the vicinity of the transport threshold, the length 
over which the number of grains transported relaxes to 
its saturated state is much shorter than the drag length 
- it decays as 1/uJ. Therefore, the overall saturation 
length is proportional to the drag length, as confirmed 

by direct measurement^^. 

The case of bed-load is still under debattPllIlIl. We 
derive here the saturation time and the saturation length 
in the simple bed-load model detailed above. As the mov- 
ing grains form a surface layer of thickness d, the number 
of moving particles per unit area adapts immediately to 
a change of shear velocity. By contrast, the grain veloc- 
ity relaxes to its asymptotic value with a characteristic 
time. This is what gives the saturation time. Neglect- 
ing the dependence of the drag coefficient on the particle 
Reynolds number, the horizontal component of the grain 
equation of motion reads: 



dvF 
~dt 
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Linearising this equation around the asymptotic value, 
we obtain the following expression for the saturation 
time: 



(51) 



Using expression (36) for Ud and typical values for the 
various parameters, we get Tg^t on the order of few y/d/g. 
The saturation length is then the length over which the 
grain moves during Tsat at velocity uP; 
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(52) 



Inserting again typical numbers in this expression, we 
get, for u close to u^, a value for Lgut on the order of few 
grain diameters. 



B. Exchange time vs saturation time 

An important problem that cannot be tackled using 
the simple transport models presented here (or any Eu- 
lerian continuous model) is the exchange between the mo- 
bile and the static phases. Such models do not aim to 
describe the Lagrangian paths of individual grains. In 
particular, recent studies have focused on the character- 
istic time a given grain spends in the transport layer be- 
fore being trapped by the becP^E^ISll, xhis time, noted 
Tox hereafter, is either called the deposition time or the 
exchange time. It is relevant in geology as it reflects 
the time scale associated with storage and reworking of 
sediments. The exchange time should not be confused 
with the saturation time. Imagine for instance the case 
where all the grains in the transport layer would move 
with a uniform and perfectly horizontal velocity. Then 
there would be no exchange with the static phase and Tox 
would be infinite, although transport could reach satu- 
ration after a very short time. Despite this conceptual 
difference, the formalism proposed irP3 leads to an iden- 
tity between the exchange time and the saturation time. 

Using our granular based transport simulations, we ad- 
dress this issue for bed load {pp/pf = 2) by tracking all 
grains with velocities above a certain value at t = 0. For 
the sake of this discussion, we have chosen this value to 
be y/gd/2, which allows us to determine the grains inside 
the transport layer at the initial time. Noting this parti- 
cle ensemble f , we define the density rit of moving grains 
at time t ~ that remains in the transport layer after a 
time t: 
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(53) 



The time evolution of nt in our numerical simulations is 
displayed in Fig. 



(50) 

an exponential re 



14 1. It can be seen that nt{t) follows 



axation with time 



nt{t) ^ {no - noo) exp {-1/%^) + n^ 



(54) 



where ng and the initial and asymptotic values, 

respectively. As the grains that no longer move have been 
exchanged with the static phase, the relaxation time of 
nt is by definition the exchange time Tox. 
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FIG. 14. (a) Time decay of the fraction of grains with an ini- 
tial velocity above that remain in the transport layer 
at time t, for different shear velocities. The solid line is the 
exponential fit (Eq. 541. Panels (b) and (c) respectively show 



the ratio of the asymptotic to initial value and the characteris- 
tic exchange time, as a function of the rescaled shear velocity. 



When analyzed for different shear velocities, the frac- 
tion of grains re-entrained in the flow after being trapped 
by the bed, which is given by the ratio noo/^^o (Fig. 14 d), 
depends weakly on 



'Q. However, it seems to tend to 
zero at the threshold, which is a reasonable limit as all 
transported grains should be eventually trapped by the 
surface and replaced by new ones. The exchange time is 
also roughly constant, with a mean value Tex — 70-y/c?/.g 
(Fig. 14 3). This time is larger at least by one order of 
magnitude than the saturation time. This means that 
exchange between the bed and the transport layer is not 
the dominant mechanism for the relaxation of the sedi- 
ment flux towards saturation. 



VI. CONCLUSIONS 

The aim of this paper was to present a novel numerical 
approach for sediment transport based on a discrete el- 



ement method (DEM) for particles coupled to a contin- 
uum Reynolds averaged description of hydrodynamics. 
We have studied the effect of the grain to fluid density 
ratio Ppl Pf and showed that we can reproduce both (sub- 
aqueous) bed load at Pp/ Pf close to unity, where trans- 
port occurs in a thin layer at the surface of the static 
bed, and (aeolian) saltation at large pp/py, where the 
transport layer is wider and more dilute. 

We have studied the mechanisms controlling steady, or 
saturated transport. In the bed load case, saturation is 
reached when the fluid borne shear stress at the interface 
between the mobile grains and the static grains is reduced 
to its threshold value. The number of grains transported 
per surface unit is therefore limited by the available mo- 
mentum at the bed surface. However, the fluid velocity 
in the transport layer remains almost undisturbed so that 
the mean grain velocity scales with the shear velocity u.^, . 
In the saltation case, particles in motion are able to eject 
others when they collide with the static bed, and satura- 
tion is reached when one grain is statistically replaced by 
exactly another one after collision. As a consequence, the 
mean grain velocity scales on the shear velocity thresh- 
old Ud, independently of u*. This provides evidence for 
a strong negative feedback of the moving grains on the 
flow within the transport layer, where the wind velocity 
is reduced. In both bed load and saltation regimes, the 
number of grains transported per unit area is found pro- 
portional to the distance to threshold O — Qdi with an 
identical prefactor on the order of 1 /cf . 

We have systematically varied the density ratio in or- 
der to reveal the transition between these two trans- 
port regimes. This is also relevant for sediment trans- 
port i n ext raterrestrial atmospheres (Mars, Venus and 
Titan ^I^ZIMIMl, have shown that the properties of bed 
load transport are observed when Pp/pf ^ 10), whereas 
those of aeolian saltation are well established when Pp/pf 
is larger than few hundreds. Finally, we have discussed 
the saturation transient of sediment transport. Based on 
the mechanisms identified in the steady case, we have 
derived expressions for the saturation time and length in 
the two regimes. In the bed load case, we have also shown 
that the exchange time, which reflects the time scale as- 
sociated to exchange of particles between the mobile and 
static phases is an order of magnitude larger than the 
saturation time. 

This study could be continued in different directions. 
First, it would be interesting to look at the case where the 
bed is non erodible. This situation has been experimen- 
tally investigated in the aeolian regime^^, showing a much 
wider transport layer A, and new scaling laws for A, the 
roughness zq and the flux (jsat as a function of . Further 
work should also be done to perform direct measurements 
of Lsat and Tsat- However, the study of inhomogeneous 
or unsteady situations requires a flner implementation of 
the model, especially for averaging procedures. Finally, 
one shall take into account the turbulent fluctuations and 
to address the case of suspended transpor1)2Si. 
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